Deep inner-shell multiphoton ionization by intense x-ray free-electron laser pulses 
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We have investigated multiphoton multiple ionization dynamics of argon and xenon atoms using 
a new x-ray free electron laser (XFEL) facility, SPring-8 Angstrom Compact free electron LAser 
(SACLA) in Japan, and identified that Xe"^ with n up to 26 are produced predominantly via 
four-photon absorption as well as Ar""*" with n up to 10 are produced via two-photon absorption 
at a photon energy of 5.5 keV. The absolute fluence of the XFEL pulse, needed for comparison 
between theory and experiment, has been determined using two-photon processes in the argon atom 
with the help of benchmark ab initio calculations. Our experimental results, in combination with 
a newly developed theoretical model for heavy atoms, demonstrate the occurrence of multiphoton 
absorption involving deep inner shells. 



Multiphoton processes in the optical regime are well- 
known phenomena investigated for decades. The ad- 
vent of extreme ultraviolet (EUV) Q and x-ray [1| 
free-electron lasers (FELs), with femtosecond pulse 
widths, has led to renewed interest in multiphoton pro- 
cesses in the EUV to x-ray spectral region. See, for 
example, some of the recent works on atoms and small 
molecules carried out at FLASH in German y [4 -161 
and at the SCSS test accelerat or p in Japan [THllj. as 
well as at LCLS Q in USA [IMj The motivation for 
these studies has been not only to reveal these pathways 
of the multiphoton multiple ionization newl y op ened by 
these light sources (see, e.g., H i, [H ESSl) 



but also to employ these processes as the basis for a new 
typ e of spectroscopy for chemical analysis (see, e.g.. 



18l . |25|). The relevance of multiphoton multiple 



ionization processes for electronic radiation damage in 
materials has also been noted l^, 2^. Electronic radi- 
ation damage due to multiphoton processes is a crucial 
issue for x-ray imaging using an XFEL. So far, how- 
ever, multiphoton experiments have been limited to the 
photon energy range up to 2 keV, i.e., the upper pho- 
ton energy limit of the atomic, molecular, and optical 
physics beam line at LCLS. 



In March 2012, a new XFEL facility, the SPring-8 
Angstrom Compact free electron LAser (SACLA) [27| . 
started user operation in Japan. Using this new facility, 
we have investigated multiphoton multiple ionization 
dynamics of argon and xenon atoms in intense hard 
x-ray pulses. The Ar if-shell thresholds are around 
3 keV, and the Xe L-shell thresholds are around 5 keV. 
The intense x-ray pulses from SACLA, with a photon 
energy above those deep inner-shell thresholds, induce 
complex multiple ionization dynamics characterized by 
the absorption of several photons. Wc identify that for 
a single Xe atom, absorption of 4 photons of 5.5 keV in- 
duces emission of up to 26 electrons. As a consequence 
of intraatomic electron-electron interactions, each pho- 
ton causes the ejection of more than 6 electrons on av- 
erage. 

From the pioneering work on the light neon atom, we 
learned that x-ray multiphoton multiple ionization is 
well characterized by a sequence of one-photon ioniza- 
tion accompanied by decay processes [lJ,i2(| . The com- 
plexity increases for heavy atoms, where one-photon 
ionization of a deep inner shell initiates a decay cas- 
cade, i.e., a series of decay steps leading to the emission 
of several electrons [i^ . The extremely large number of 
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FIG. 1: (Color online) Experimental configuration. 



x-ray photons in an XFEL pulse is responsible for trig- 
gering further photoionization from deep inner shells 
during and after a decay cascade. Thus sequential mul- 
tiphoton absorption involving deep inner shells becomes 
very complicated. As a consequence of its sequential 
nature, x-ray multiphoton absorption depends primar- 
ily on the fluence, which is the number of photons (or 
the pulse energy) per unit area, rather than the (peak) 
intensity [itI Ii9| . The minimum fluence required for 
multiphoton absorption is estimated as follows. For 
example, the photoionization cross section of neutral 
Xe at 5.5 keV is ^0.166 Mb, so the fluence to satu- 
rate one-photon absorption is ^6x10^° photons/yum^ , 
which corresponds to ^50 /iJ//im^ for 5.5-keV photons. 
If the fluence of an x-ray pulse is close to or higher than 
this fluence, one expects that more than one photon will 
be absorbed by a xenon atom within one pulse. 

The experiment has been carried out at the exper- 
imental hutch 3 (EH3) of the beam line 3 (BL3) of 
the SACLA in Japan The photon energy was 

set at 5.5 kcV. The photon band width was '-^60 eV 
(FWHM). The repetition rate of the XFEL pulses was 
10 Hz. The pulse width has not been measured, but 
was estimated to be in the range between 10 and 30 fs 
(FWHM). Shot- by-shot pulse energy was measured by 
the beam-position monitor [29| located upstream of the 
beam line described below. The monitor was calibrated 
by a calorimeter so that output signals from the 
monitor could be transformed to the absolute value of 
the pulse energy. The measured values were 239 fj.J per 
pulse on average. 

Figure [T] shows the experimental configuration. A 
Kirkpatrick-Baez (KB) mirror system is permanently 
installed at EH3 [31| . The focal length (work distance) 
is ~1.3 m. The XFEL beam is focused by the KB 
mirror system to a focal size of ~1 /xm (FWHM) in 
diameter. The Rayleigh length is ~8 mm. The sample 
gas (argon or xenon) was introduced as a pulsed super- 
sonic gas jet [sll to the focus point of the XFEL pulses, 
in the ultrahigh-vacuum reaction chamber. The gas 
beam at the reaction point was estimated to be ~2 mm 
(FWHM) in diameter. Thus, the source volume of the 
ions was roughly a cylindrical shape of ^1 fim in diam- 
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FIG. 2: The ion time of flight spectrum of Xe recorded at 
the photon energy of 5.5 keV at full XFEL pulse energy. 



eter and ~2 mm along the XFEL beam. After crossing 
the gas jet at right angles, the XFEL beam came out 
of the vacuum chamber via a beryllium window. The 
relative x-ray pulse energy was measured shot-by-shot 
by a p-intrinsic-n (PIN) photodiode, after the pulse en- 
ergy was reduced by aluminum sheets of 0.2-0.5 mm 
thickness so that linear response of the photodiode was 
assured. The shot-by-shot pulse energy fluctuation was 
±25% (50% FWHM). 

Ions produced in the source volume described above 
were extracted towards the ion time-of-fiight (TOF) 
spectrometer [33| equipped with microchannel plates 
(MCP) and a delay-hne anode (Roentdek HEX80) [si]. 
Signals from the MCP and the delay-line anode were 
fed to an 8-channel digitizer. The wave forms recorded 
by the dig itizer were analyzed by a software discrim- 
inator 



35| and the arrival time and the arrival posi- 



tion of each ion were determined. The voltages on the 
spectrometer are tuned for the best mass resolution. 
The simulation of ion trajectories shows that the TOF 
slightly depends on the ion arrival position on the de- 
tector, as a consequence of different departure positions 
perpendicular to spectrometer axis. A compensation on 
the measured TOF is introduced to improve the mass 
resolution further. Figure [2] depicts the TOF spectrum 
for xenon ions after compensation of the TOF. We can 
clearly see ions with a charge state up to -1-26 and well- 
resolved isotopes at different charge states. 

Let us now describe our theoretical approach to x- 
ray multiphoton multiple ionization dynamics. We em- 
ploy the XATOM toolkit , based on the rate equation 
approach and the Hartree-Fock-Slater method. 
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For Ar, we use the approach described m Ref. 23|, 
adding shakeoff processes within the sudden approx- 
imation [3^. For Xe, even this straightforward se- 
quential ionization model becomes tremendously chal- 
lenging, because (a) there are too many electronic 
configurations involving multiple holes and too many 
atomic data (photoionization cross sections, Auger and 
Coster-Kronig rates, and fluorescence rates), and (b) 
the matrix size for the set of coupled rate equations 
is too large to be solved directly. We have addressed 
the latter by introducing a Monte Carlo approach in 
Refs. [l^lij] where for Xe M-shell ionization dynamics 
~10^ coupled rate equations were solved and ^4x10'' 
atomic data were pre-calculated. For Xe L-shell ioniza- 
tion dynamics considered here, however, the complexity 
is further increased by a factor of 21. thus it becomes 
formidable to pre-calculatc all the atomic data required 
for the rate equations. Therefore, we extend xatom by 
applying the Monte Carlo procedure for both calculat- 
ing atomic data and searching probable ionization path- 
ways. When a certain configuration is visited during 
an ionization pathway chosen by the Monte Carlo sam- 
pling, atomic data are calculated for the corresponding 
configuration. In this way, atomic data are computed 
only when they are required. This extension enables 
us to treat ionization dynamics of heavy atoms with no 
limit of configurational space. Saving in computational 
time is dramatic. For example, one calculation of Xe at 
5.5 keV and 190 /iJ//im^ takes 2 days using one CPU 
on the laboratory workstation, which would take more 
than 3 years if all atomic data had to be calculated. 

In the present calculations, the photon energy is fixed 
at 5.5 keV. The pulse shape is assumed to be a Gaus- 
sian of 30 fs (FWHM). In the regime of sequential mul- 
tiphoton absorption, the results are not sensitive to the 
pulse shape or spikiness of individual pulses [l^]. We 
also did the calculation for 10 fs (FWHM) at a sin- 
gle XFEL fluence and confirmed that the results are in 
good agreement. All calculated results presented here 
are obtained by 3-dimensional integration [3] over the 
interaction volume of 1 //m x 1 ^m x 2 mm (all lengths 
are given in FWHM), according to the instrumental 
configuration. 

In order to compare the theoretical results with the 
experimental results, we need the peak fiuence of the 
XFEL pulse employed for our experiment. The fluence 
at the center of the x-ray beam focus defines the peak 
fluence i^poak, which is given by 

4 In 2 E 

-Fpcak = X — X T, 

where E is the nominal pulse energy from the monitor, 
T is the transniission(%), A is the focal area given in 
FWHM X FWHM, and the coefficient of 41n2/7r comes 
from the assumed Gaussian focal shape. In order to es- 
timate the peak fiuence of the XFEL pulse, we have 
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FIG. 3: (Color online) Experimental and theoretical charge 
state distributions of Ar at the photon energy of 5.5 keV. 



measured the charge state distribution of Ar atoms 
and compared it with our benchmark calculations. At 
5.5 keV, the /-iT-shell electrons of Ar can be ionized. The 
number of all possible multiple-hole configurations, or 
equivalently the number of coupled rate equations is 
1323, which is easily solvable without the Monte Carlo 
method. 

The experimental and theoretical charge state distri- 
butions of Ar are shown in Fig. [3] Theoretical charge 
state distributions are scaled so that the sum of the 
individual charge state yields is equal to the total ion 
yield in the experiment. We find that the theoretical 
distributions up to the charge state of 4-6 do not vary 
when varying the XFEL fluence, illustrating that these 
ions arc produced by single photon absorption. In gen- 
eral, they are in good agreement with experiment, ex- 
cept that in theory the Ar'^+ is slightly overestimated 
while Ar^+ and Ar'''*' are underestimated. Further, 
we flnd that both experimental and theoretical yields 
of Ar^+ and Ar^+ depend quadratically on the XFEL 
fluence. We thus can estimate relative contributions 
from two- and one-photon processes using the ratio of 
R = (Ys++Y9+)/(Y3++Yi+), where r„+ is the yield of 
Ar"+ . Fitting the theoretical value of R to the exper- 
imental one by the least squares method, we have ob- 
tained 70 /^J//im^ for the peak fiuence. The uncertainty 
of this estimate mainly stems from our benchmark the- 
oretical calculations. We expect that this uncertainty 
is ^±10%. The resulting value of 70±7 fii/ ^ir? is con- 
sistent with the empirical knowledge from other exper- 
iments, that the transmission of the KB mirror system 
is ^50% and that the focal area is ^2 fim^, though 
these numbers may include uncertainties of a factor of 
two. 

Figure |4] depicts the charge state distribution of Xe 
at the peak fluence of 70 /iJ/^m^, which is determined 
via the calibration using Ar. The charge state distribu- 
tion varies as the peak fluence varies (not shown here) . 
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FIG. 4: (Color online) Experimental and theoretical charge 
state distributions of Xe at the photon energy of 5.5 keV 
with the peak fluence of 70 fiJ/nm^ . 



We also compare the theoretical charge state distribu- 
tions at this peak fluence with experiment. The pho- 
ton energy of 5.5 keV is above the L-shell threshold for 
charge states up to -1-23 according to our calculations. 
There is no signature of resonance-enabled ionization 
enhancement [19[ because the fluence is not enough to 
form high charge states that initiate resonance excita- 
tion. The discrepancy between theory and experiment 
may be attributed to the non-relativistic treatment and 
lack of shakeoff in the current theoretical model. The 
shakeoff process can further ionize valence electrons af- 
ter photoionization, and some decay channels might be 
absent without relativity Inclusion of both rel- 

ativity and shakeoff tends to produce higher charge 
states. Thus, in the current model, the formation of 
high charge states is somewhat suppressed in compar- 
ison with the experimental results. In spite of these 
limitations of the current model, the experimental and 
theoretical results are in reasonable agreement, at least 
semiquantitatively. 

To obtain the fluence dependence of the yields for the 
xenon ions at different charge states in a pulse energy 
range wider than the shot-by-shot pulse energy fluctua- 
tion, we recorded the spectra not only at the full XFEL 
pulse energy but also at a pulse energy attenuated by 
an aluminum foil of 25 fim thickness located upstream 
of the KB mirror system (see Fig. [1]). The attenuated 
pulse energy was 38% of the full pulse energy according 
to the reading of the PIN photodiode. In Fig. [5l the 
ion yields for Xe""*" (n=8, 14, 18, and 24) are plotted as 
a function of the peak fluence. To obtain several data 
points with respect to the peak fluence, we first merged 
the results measured at two different pulse energies, and 
then re-binned the data. Straight lines with a slope of 1, 
2, 3, and 4 are also shown as a guide to the eye. The flu- 
ence given on the horizontal axis is on the absolute scale 
determined by the Ar calibration. All theoretical ion 
yields are scaled by a single factor such that the yield of 
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FIG. 5: (Color online) XFEL fluence dependence of the ion 
yields for Xe"+ (n=8, 14, 18, and 24). Closed circles with 
the error bars depict the experimental results and sohd lines 
depict the theoretical results. Lines with slope p = 1, 2, 3, 
and 4 are also shown as broken Hnes to guide the eye. The 
uncertainty of the peak fluence is expected to be ~±10%. 



Xe^~^ at 70 /iJ /^m^ is matched with that in the exper- 
iment. The yield of Xe*"*" exhibits a slope of less than 
one, illustrating that this ion is produced by single- 
photon absorption and that the single-photon absorp- 
tion process is close to saturation. The yield of Xe^^"*" 
exhibits a slope close to two, the yield of Xe^®"*" exhibits 
a slope close to three, and the yield of Xe^'*+ exhibits a 
slope close to four. These slopes directly suggest that 
these ions are predominantly produced by two-, three- 
, and four-photon absorption, respectively. Although 
theory reproduces the charge state distributions only 
semiquantitatively, it reproduces very well the fluence 
dependence of the individual charge states. This agree- 
ment confirms the number of absorbed photons yielding 
the individual charge states and thus fully confirms that 
L-shell multiphoton processes of the xenon atom take 
place in the x-ray regime at 5.5 keV. 

It is worth noting the relevance of the present work 
to the emerging area of femtosecond crystallography 
with XFELs [33, [s^- Electronic radiation damage, 
especially to heavy atoms, is inevitable during fem- 
tosecond XFEL pulses, as seen in the present work 
on the xenon atom. The high-inten sity version of 
multiwavelength anomalous diffraction |26| beneficially 
exploits multiphoton multiple ionization dynamics of 
heavy atoms embedded into macromolecules and pro- 
vides a new path to determine macromolecular struc- 
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ture using XFELs. Therefore, a comprehensive picture 
of multiphoton multiple ionization dynamics is crucial 
to success in femtosecond x-ray imaging with XFELs. 

In summary, we have investigated multiphoton mul- 
tiple ionization dynamics of the xenon atom using the 
new XFEL facility SACLA as well as a newly devel- 
oped theoretical model, and identified that Xe"'^ with 
n > 24 arc produced predominantly via four-photon ab- 
sorption, demonstrating the occurrence of multiphoton 
absorption of a single atom in the x-ray regime above 
2 kcV for the first time. The absolute fluence of the 
XFEL pulses has been determined using two-photon 
processes in argon atoms with the help of benchmark 
ab initio calculations. 

The experiments were performed at SACLA with the 
approval of JASRI and the program review committee 
(No. 2012A8036). This study was supported by the 
X-ray Free Electron Laser Utilization Research Project 
and of the Ministry of Education, Culture, Sports, Sci- 
ence and Technology of Japan (MEXT), by the Japan 
Society for the Promotion of Science (JSPS), and by 
the IMRAM project. 
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